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Abstract 

We develop further the proposal of arXiv:0806.3692 that a new state of 
matter - charged condensate of spin-0 nuclei - may exist in helium-core dwarf 
stars. The charged condensate and its fluctuations are described by an effec- 
tive field theory Lagrangian. The spectrum of bosonic fluctuations is gapped, 
while electrons, at temperatures of interest, give rise to gapless excitations 
near the Fermi surface. These properties determine the evolution of the dwarfs 
with condensed cores. In particular, we show that such dwarf stars would cool 
significantly faster than their crystallized counterparts. As a result, the lu- 
minosity function for the helium-core dwarfs will have a sharp drop-off after 
the condensation. It is tempting to interpret the recently discovered abrupt 
termination of a sequence of 24 helium-core dwarf candidates in NGC 6397 
gnature of the charged condensation. 



1. Introduction and Summary 



There is roughly five orders of magnitude between what may be called the atomic 
scale, Angstrom ~ 10~^ cm, and the nuclear scale, Fermi ~ 10~^^ cm. If in a neutral 
system of a large number of electrons and nuclei average inter-particlc separations 
are between these two scales, then neither atomic nor nuclear effects will play any 
significant role. Moreover, the nuclei can be treated as point-like particles. 

Such a system of electrons and nuclei constitutes cores of white dwarf stars. Up 
to a factor of a few, these are roughly Earth's size solar-mass objects; their central 
mass density may range over ~ (10^ — 10^'') g/cm?, most of them being near the 
lower edge of this interval. Since the dwarf stars exhausted thermonuclear fuel in 
their cores already, they evolve by cooling [1]; the ones that we consider in this work 
cool from ~ 10^ K down to lower temperatures. 

At temperatures of interest the electron de Broglie wavelengths overlap with 
each other, and they obey the laws of quantum statistics. Hence, the electron 
properties are universal throughout all the dwarf stars; their number density may 
range in the interval ~ (0.1 — 5 MeVY, and at temperatures of interest they form 
a degenerate quantum gas (their Fermi energy is greater that their interparticle 
Coulomb interaction energy). 

As to the nuclei - since they're heavier and their de Broglie wavelengths are 
shorter - they may exhibit different properties depending on the values of density 
and temperature at hand: they may crystallize [2] as it should be the case for 
a majority of dwarf stars, or alternatively, may form a quantum liquid [3], when 
densities are high enough. Identification and qualitative description of quantum 
liquid signatures in the cores of high-density cool white dwarf stars is a main subject 
of the present work. 

As a typical dwarf star cools down, the Coulomb interaction energy in a classical 
plasma of charged nuclei will significantly exceed their classical thermal energy, and 
the nuclei, in order to minimize energy, would organize themselves into a crystal 
lattice [2] . In most of these cases quantum effects of the nuclei should be negligible; 
for instance, the Debye temperature should be less than the temperature at which 
crystallization takes place, and the de Broglie wavelengths of the nuclei should be 
much smaller than the average intcrnuclcar separations. This indeed is the case in 
majority of white dwarf stars, the cores of which are composed of carbon and/or 
oxygen nuclei and span the interval of mass densities around ~ (10^ — 10^) g/crri^. 

However, there exists a class of dwarf stars in which the nuclei enter the quan- 
tum regime before the above- described classical crystallization process sets in [4, 5]. 
Among these, furthermore, there is a relatively small subclass of the dwarf stars for 
which the temperature Tg, at which the de Broglie wavelengths of the nuclei begin 
to overlap, is higher than the would-be crystallization temperature Tcrj/st (see [3] 
and discussions in Section 2 below). Then, right below Tc, the quantum-mechanical 
uncertainty in the position of the charged nuclei is greater that the average inter- 
nuclear separation. This is diametrically opposite to the crystallized state where the 
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nuclei would have well- localized positions with slight quantum-mechanical fuzziness 
due to their zero-point oscillations. 

What is then an adequate description of such a state? It was argued in [3, 6] 
that such a system, instead of forming a crystaUine lattice at Tcryst, could condense 
at Tc > Tcryst, owing to the quantum-mechanical probabilistic "attraction" of Bose 
particles to occupy one and the same zero-momentum state, and leading to a quan- 
tum liquid in which the charged spin-0 nuclei would form a macroscopic quantum 
state with a large occupation number - the charged condensate. 

It was proposed in [3] that the conditions in certain high-density helium-core 
white dwarfs (He WDs) arc appropriate to form such charged condensate. Here we 
will argue that similar effects could take place in lower-density He WDs, which are 
more relevant for observations [7, 8, 9], as well as in superdense (~ 10^° g/crn?) 
carbon-core white dwarfs (C WDs). 

The observational consequences of the charged condensate in astrophysical ob- 
jects can be significant. The bosonic part of the liquid is superconducting as its 
spectrum of long- wavelength fluctuations exhibits a mas gap; hence these quasi- 
particles contribute to the specific heat of the substance only in an exponentially 
suppressed way [3]. This affects the cooling rate for dwarf stars. In particular, we 
will show that the He WDs with the condensed cores would cool much faster than 
the ones in which condensation does not take place. 

We should emphasize though that a careful consideration of white dwarf cool- 
ing rates requires detailed account of composition and dynamics of their envelopes, 
which is a complex and less certainly known subject (for detailed discussions see, 
e.g., [10, 11]). We will not attempt here to enter these studies. Instead, we con- 
sider an over-simplified model just to emphasize our main point that if the charged 
condensation takes place in the cores of white dwarfs, they'll necessarily cool faster, 
and there will be a drop-off in their luminosity function. 

One may wonder whether the effects that we're discussing can be entirely obfus- 
cated by the uncertainties in the envelope composition. Although this may well be 
the case for cooling times, there is nevertheless a model independent prediction of 
the charged condensation: the luminosity function will have a sharp drop-off after 
the condensation, with subsequent growth governed by a shallower slope. 

The recently discovered and studied in Ref. [9] 24 hehum-core dwarf candidates 
in NGC 6397 exhibit a termination of the sequence at low luminosities. It is tempting 
to attribute this to the drop-off due to the charged condensation. Whether this 
proposal can withstand more detailed scrutiny remains to be seen. 

The organization of the paper is as follows: In Section 2 we give qualitative 
arguments why condensation may be preferred over crystallization in certain WDs. 
In Section 3 we summarize and further develop the effective field theory approach 
of [3, 12] to charged condensation, and describe the spectrum that is relevant for 
cooling of WDs. In Section 4 we calculate cooling times for He WDs and compare 
them with those of core-crystallized WDs. Last but not least, we discuss the drop-off 
of the luminosity function after charged condensation. 
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2. Condensation versus Crystallization 



We concentrate on WDs that below certain temperature cool by releasing stored 
in them heat from their surface. At the beginning of that stage the ions (nuclei) 
form a classical Bose gas that cools with constant specific heat (the so-called Mestel 
cooling [1]) down to some temperature Tcryst, at which Coulomb repulsion starts to 
dominate by about two orders of magnitude over the average thermal energy, and 
the core undergoes crystallization transition [2]. The crystallization temperature in 
the classical regime is quantified by the ratio 

^Coulomb _{Zey/4:TTd 
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where e is the electric charge, Ze is the charge of a nucleus, denotes the Boltz- 
mann constant, Jq is the electron number-density, d = (SZ/AnJoY^^ is the average 
inter-ion separation, while T denotes the core temperature of the star. The ion 
plasma is expected to undergo crystallization once the temperature drops to a value 
for which T ~ 180 [13, 14, 15]. 

The above arguments are entirely classical. The temperature scale that deter- 
mines the classical versus quantum nature of the crystallization transition is the 
Debye temperature 



e,^^, (2) 
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where Qp is the plasma frequency of the ion gas, Jh — Jq/Z denotes the ion num- 
ber density and mn denotes the mass of a single ion (the subscript H stands for 
"heavy"). Up to a factor of ^/S, Hp coincides with the frequency of zero-point 
oscillations of the ions in crystal sites, wq = Vtp/ v^. 

The white dwarfs with Tcryst > do cool according to the above- described classical 
scenario. Often however, 6d may significantly exceed Tcryst [4]. In such a case, 
quantum zero-point oscillations should be taken into account in order to derive the 
crystallization temperature. This seems to delay the formation of quantum crystal, 
lowering Tcryst from its classical value at most by about ~ 10% [5]. Since this is a 
small change for the estimates that we're after here, we will consider the classical 
value of Tcryst to be a good approximation, even in the quantum case, keeping in 
mind that Tcryst may overestimate somewhat the crystallization temperature of the 
substance. 

However, there exists a third and very important temperature scale, relevant for 
studying the cooling of white dwarf interiors. It is the "critical" temperature Tc, at 
which the de Broglie wavelengths of the ions start to overlap 

47r2 
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Below Tc quantum-mechanical uncertainties in the ion positions become greater than 
an average inter-ion separation. Hence the latter concept looses its meaning as a 
microscopic characteristic of the system, and the ions enter a quantum-mechanical 
regime of indistinguishability^ . Below Tc the wavefunction of the many-body system 
of spin-0 ions should be symmetrized, and this would unavoidably lead to proba- 
bilistic "attraction" of the bosons to occupy the same quantum state. 

Therefore, when crystallization temperature Tcryst is lower than Tg, the system 
may instead undergo condensation into a macroscopic zero-momentum quantum 
state with large occupation number - the charged condensate. 

Once in the condensate, the boson positions are entirely uncertain while their 
momenta equal to zero. In order for such a system to crystallize later on, each of 
the bosons should acquire the momentum determined by the zero-point energy of 
the crystal ions, kg ~ 2mH^p/ "n/S- The latter is greater than the Fermi momentum, 
as well as typical momenta of fermionic quasiparticles. Hence, the fermions will 
not be able to transmit to the condensed bosons momenta comparable with |ko|. 
Therefore, the transition to the crystallized state can only happen spontaneously. 
Such a transition could take place as it would lower the energy of the entire system 
due to the favorable electrostatic screening. However, given that the spectrum of 
bosonic quasiparticles is gapped, the process will be exponentially suppressed at 
temperatures below the gap scale. Hence, the condensate should not be expected 
to undergo subsequent crystallization, at least for a long period of time. 

Another crucial question is whether the expression for the critical temperature 
(3) gives an accurate estimate for the actual condensation temperature Tcond at which 
the phase transition would take place. If we were to deal with a non-interacting sys- 
tem of Bosons then the known BE condensation temperature, Tj,^j^ — {1.27/mH(P), 
would have been an order of magnitude smaller than what the estimate (3) sug- 
gests. However, it has been known that already weak repulsive interactions between 
bosons increase the condensation temperature; this is consistent with ones expecta- 
tion that the repulsion makes easier for the condensation in the momentum space to 
take place (indeed, the BE condensation is a condensation in the momentum space, 
while the coordinate space wave- functions are entirely delocalized) see, e.g., [16] and 
references therein. 

In the case of weakly interacting bosons the increase of the condensation tem- 
perature is small since the interactions are weak. In our case, however, interactions 
between spin-0 nuclei are strong, in a sense of a many-body system, as we are about 
to argue below. In this case we would expect Tcond ^ ^cwkI- Since we have no means 
to evaluate Tcond accurately, we use the expression (3) as a reasonable and physically 
motivated approximation for the interacting system Tc ~ Tcond- To this end, the 

^The de Broglie wavelength above is defined as A^b = 27r/|k|, where k^/2m// = 2>kBT/2. We 
define as the temperature at which X^b — d. Note that this differs by a numerical factor from 
the standard definition of the therniMl do Broglie wavelength. A = ^J2'K /mkBT, that appears as a 
natural scale in the partition function. See comments on the rationale for our choice of below. 
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criterium that we adopt for the condensation to take place is 

Tc {a few) nrysf (4) 

The charged condensate is somewhat similar, but also differs by its strong cou- 
pling, from the Bose-Einstein condensate of charged spin-0 particles. Analytic stud- 
ies of a BE condensate of charged scalars in the past (see, e.g., Ref. [17] and refer- 
ences citing it) relied on a small departure from the condensation of free particles, 
and used an expansion in a parameter that is defined as the ratio of the average 
intcrparticle separation to the would-bc Bohr radius for the boson. For a weakly 
non-ideal system of bosons the rs-ratio is small <^ 1, and the ground state is 
the weakly-non-ideal BE condensate. Furthermore, as zero-temperature numerical 
simulations show (see, e.g., [18, 19], and references therein), for 1 ^ ^ 160 the 
condensate and crystal state begin to coexist. Furthermore, for > 160 the entirely 
crystallized state is a ground state. Presumably, the finite temperature effects would 
increase the value of at which the crystallization takes place only by a factor of a 
few, but not by an order of magnitude (see, e.g., [20]). 

In the system of nuclei and electrons the r^-ratio cannot be made small without 
entering the regime where nuclear interactions become dominant. Instead, in the 
case considered in [3], ~ (10 — 100). In the present work we will consider even 
larger values, ~ 10^, which are relevant for observed He WDs. Yet we argue 
that because the temperature Tc at which the de Broghe wavelengths of bosons 
start to overlap is greater than the crystallization temperature, the system, upon 
gradual cooling, should settle in the charged condensate state due to the quantum- 
mechanical probabilistic properties of the indistinguishable bosons. In fact, our 
criterium > T^ryst implies that the charged condensation could take place at fine 
temperature Tc for Ts < 2400. It is also worth pointing out in this regard that 
the expansion parameter in our case ends up being l/r^, as long as the system is 
described in an effective Lagrangian approach. 

The seeming contradiction with the numerical results may be reconciled by the 
fact that the charged condensate can only be a metastable state [21], while the crystal 
should arguably be the lowest energy state for » 160. Most of the numerical 
simulations use the test wavefunction approach that minimizes the energy; finding 
a long-lived state which represents only a local minimum may not be easy in this 
approach^. 

One check of this proposal is that the small fiuctuations of the charged condensate 
have no unstable modes, suggesting that it represents at least a local minimum. Also, 
these fluctuations are rather different from those of a would-be crystal, ordinary 
clod plasma, or weakly-coupled BE condensate. For discussions of non-perturbative 
(meta) stability of the charged condensate in a different context, see [21]. 



^It is interesting to note that in Ref. [18] a metastable quantum liquid branch was observed 
for Vg > 160. At this stage it is hard to speculate whether this is the branch that we're discussing 
here. 
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3. Charged Condensate and its Fluctuations 



We use an effective Lagrangian description to study the charged condensate^. 
We focus on the zero-temperature hmit, even though reahstic temperatures in He 
WDs are well above zero (for calculations of the finite temperature effects, see [22]). 
The validity of the zero-temperature approximation is justified a posteriori and goes 
as follows: the spin-0 nuclei undergo the condensation to the zero-momentum state; 
while they do so they cannot excite their own phonons since the latter are gapped 
with the magnitude of the gap being greater than the condensation temperature. 
On the other hand, the condensing charged bosons can and will excite thermal 
fluctuations in the fermionic sector that is gap-less. Therefore, all the thermal 
fluctuations will end up being stored in the fermionic quasiparticles near the Fermi 
surface. For the later, however, the finite temperature effects aren't significant since 
their Fermi energy is so much higher, T/Jq^^ -C 10~^. 

After these comments we turn to the effective Lagrangian. By mn we denote the 
mass of a (heavy) nucleus of charge Ze and atomic number A (helium-4, carbon, 
oxygen), by /x/ the electron chemical potential, and by rUe the electron mass. The 
following hierarchy of scales, mn ^ max[/i/, me], is a starting point for the effective 
Lagrangian construction. 

We begin at scales that are well below the heavy mass scale tuh, but somewhat 
above the scale set by max[/^j, mg]. Hence the electrons are described by their Dirac 
Lagrangian, while for the description of the nuclei we will use a charged scalar order 
parameter ^{x). As it was shown in [12], in a non-relativistic approximation for the 
nuclei, an effective Lagrangian proposed by Greiter, Wilczek and Witten (GWW) 
[23] in a context of superconductivity, is also applicable for the description of the 
charged condensation, given that an appropriate reinterpret ation of its variables and 
parameters is made. 

The construction of the GWW effective Lagrangian is based on the following 
fundamental principles: it is consistent with the translational, rotational, Galilean 
and the global U{1) symmetries, preserves the algebraic relation between the charged 
current density and momentum density, gives the Schrodinger equation for the order 
parameter in the lowest order, and is gauge invariant [23]. Combined with the 
electron dynamics the GWW effective Lagrangian reads (we omit for simplicity the 
standard Maxwell term): 

^eff = V Q(rDo$ - {D,^m - + n^ii^Di - mf)^, (5) 

where we use the standard notations for covariant derivatives with the appropriate 
charge assignments: Dq = {do — iZeAg), Dj = {dj — iZeAj), Dj^^ — dn + ieA^, while 

^An effective field theory that describes the charged condensate was discussed in Refs. [3, 6]. 
In this section we briefly describe and expand some of the results of these works. 
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V{x) stands for a general polynomial function of its argument ^. 

The coefficients of this polynomial, V{x) = '^'^=QCn{x"'/ A*"-), are dimensionful 
numbers that are inversely proportional to powers of a short-distance cutoff of the 
effective field theory. 

Once the basic Lagrangian is fixed, we introduce the electron chemical potential 
term fifip^tp , where /ij = e/.' = [(Stt^Jq)^/^ + m^]^/^. This gives a nonzero electron 
number density Jo which is related to the Fermi momentum kp = (Stt^Jq)^^^. This 
is also the only term that at the tree level sets a frame in which the electron total 
momentum is zero. The quantum loop corrections due to this term will generate 
additional Lorentz- violating terms in the bosonic sector of the theory [3, 22, 12]. 

There exists a solution of the equations of motion that follow from the effective 
Lagrangian (5). This solution takes the form [3]: 

Z|$|2 = Jo, A^ = 0, P'(0) = 1. (6) 

(We use the unitary gauge $ = |$|). The condition V'{0) = 1 is satisfied by any 
polynomial functions V{x) for which the first coefficient is normalized to unity 

V{x)^x + C2X^ + .... (7) 

The above solution describes a neutral system of negatively charged electrons of 
charge density — eJo, and positively charged scalar (helium-4 nuclei) condensate of 
charge density Ze$"'"$ = eJo- This describes the condensate and not a standard 
crystal in a long wavelength approximation where lattice inhomogeneities can be 
neglected, or ordinary cold plasma. This becomes more clear after one calculates 
the spectrum of small perturbations about the homogeneous solution and finds that 
it is rather different from the spectra of nearly free BE condensate, crystal lattice 
vibrations, or from plasma fluctuations. 

Calculation of the spectrum of small perturbations is straightforward (here we 
follow conventions of [12]). There are two transverse polarizations of a massive 
photon which propagate with the conventional massive dispersion relation 

u'^k' + m'^, m'^^. (8) 



"^In a more complete treatment one should also add to the Lagrangian terms /xjvii^*^, 
A(<i>*$)^/TO|j, and the higher dimensional operators that are consistent with all the symmetries 
and conditions that lead to (5). Here fiNn denotes a non-relativistic chemical potential for the 
scalars. These terms will not be important for the low-temperature spectrum of small perturba- 
tions we're interested in, as long as A ^ 1 and Jo <C m%. However, near the phase transition point 
it is the sign of ijnr that would distinguish between the broken and symmetric phases, so these 
terms should be included for the discussion of the symmetry restoration. We also note that the 
scalar part of (5) is somewhat similar to the Ginzburg-Landau (GL) Lagrangian for superconduc- 
tivity. However, there is a significant difference between the two. The coherence length in the GL 
theory is many orders of magnitude greater than the average interelectron separation, while in the 
present case, the "size of the scalar" $ is smaller that the average interparticle distance. 
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Moreover, there is one longitudinal mode (phonon) with the following unconven- 
tional dispersion relation 



2 2 



1 + 



^ m; + — ^ . (9) 



Here M = ^mHiri^, ml = + 5rn?, where 8w? oc e^Jo/Ep is the Debye mass 
squared. As we see, u'^dkl = 0) = m^, and the bosonic collective excitations are 
gapped. Hence, the bosonic part of the whole system represents a superconducting 
component, while the fermions, at temperatures of interest, can be regarded as the 
"normal" (dissipative) component of the quantum liquid (see more below) . 

We point out that in the approximation used above there does not appear to be a 
(9(k^) term on the r.h.s. of the dispersion relation (9). Such a term emerges (see [12]) 
in the l/m% order of the heavy mass expansion, and is proportional to — m^k^/m|^. 
This term is responsible for the "roton-like" behavior of small fluctuations discussed 
in [12]. The dispersion relation (9) is applicable for k'^ '^rn'^, which corresponds to 
distances shorter than 1/m^ - the scale that encompasses a huge number of particles. 

Note also that the phonon dispersion relation (9) is different from that of a 
conventional phonon in a crystal for which uj oc |k|, near the origin. Hence, the 
substance that we're discussing is different from a crystal, and it is also different 
from classical plasma. The following hue of arguments also suggests that the charged 
condensate differs from the ordinary BE condensate of nearly-free bosons too: In 
the limit of switched off interactions, e — » 0, the dispersion relation (9) reduces 
to a; = k^/2mff. This is nothing but a dispersion relation for a lowest excitation 
in a BE condensate of free bosons of mass ttih- However, for realistic values of 
the parameters the second (momentum dependent) term in (9) is sub-dominant, 
suggesting that the charged condensate has a significant departure from the BE 
condensate of nearly-free bosons. 

What is important for the present work is that the bosonic collective excitations 
give rise to exponentially suppressed contributions to the value of specific heat of this 
substance. Typical suppression scales as exp(— m^/T), and since for the He WDs 

~ 3 KeV, at temperatures below 10^ K, these contributions are proportional to 
~ exp(— 30) and can be neglected. 

This is not the case in dwarf starts in which the nuclei formed a crystalline 
structure. There, the dominant contribution to the specific heat comes from a crystal 
phonon. The latter has a linear dispersion relation a; oc |k|, and its contribution to 
the specific heat scales with temperature as T^. 

Therefore, to a good approximation, there is no phonon contribution in the 
charged condensate case, while it is present when the cores crystallize. 

As to the electrons, their behavior is similar in both crystal and condensate 
cases. At temperatures of interest they form a degenerate Fermi gas with gap-less 
excitations near the Fermi surface. Their contribution to the specific heat scales 
linearly with temperature. In the case of crystallized cores this is sub-dominant to 
the specific heat due to the crystal phonon. For the charged condensate, on the 
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other hand, the fcrmions are the dominant contributors to the specific heat. These 
properties make a diff^erence in coohng of dwarf starts, as we'll discuss in detail in 
the next section^. 

Finally, let us discuss briefiy the question of impurity (hydrogen, helium-3, etc.) 
nuclei that may be present in the cores of white dwarfs. The static potential between 
two impurity nuclei of charge Qi and Q2 consists of two parts [12]: 

Vstat = aemQiQ2 cos(Mr H \ . (10) 

\ r TT M®r^ / 

The first, exponentially suppressed term modulated by a periodic function, is due 
to cancellation between the screened Coulomb potential and that of a phonon. The 
second term, which exhibits a power-like behavior modulated by a periodic function, 
is due to the existence of gap-less excitations near the Fermi surface. This gives a 
generalization of the Friedel potential (see, e.g., [24]) to the case when in addition 
to the fcrmionic excitations the collective modes associated with the charged con- 
densate are also taken into account^. As we see, the potential is not sign-definite. It 
can give rise to attraction between like charges; this attraction is due to collective 
excitations of both fermionic and bosonic degrees of freedom. This represents a gen- 
eralization of the Kohn-Luttingcr [26] effect to the case where on top of the fermionic 
excitations the collective modes of the charged condensate are also contributing (for 
the discussion of associated superconductivity, see footnote 5). 

For the physical conditions present in dwarf stars, the second term in (10) is 
dominating. This term is strongly suppressed because of the phonon subtraction 
[12]. Thus, attraction between like charges due to this potential could give Cooper 
pairing between impurities (even if their concentration was significant enough) only 
at extremely small temperatures (see footnote 5). 



^In the charged condensate Cooper pairs of electrons can also be formed, however, the cor- 
responding transition temperature, and the magnitude of the gap, are suppressed by a factor 
exp(— l/e^jj), where Cg^j is proportional to the value of the inter-electron potential that contains 
both screened Coulomb and phonon exchange. The fact that this potential has attractive domain, 

but is very small, is suggested by the static potential foimd in [12] (sec eq. (10) below); the latter 
is down by a power of a large scale M. In other words, the static potential and the zero-zero 
component of the propagator are both suppressed as Dqo ~ l/vg ~ (10~^ — 10~^), where rg was 
discusscid in Sc;c;tion 2. Moreover, taking into account the frequency dependence via the Eliash- 
berg equation does not seem to change qualitatively the conclusion on a strong suppression of the 
pairing temperature. 

Hence, even though the bosonic sector (condensed nuclei) is superconducting at reasonably high 
temperatures <, 10^ K, interactions with gap-less fermions could dissipate the superconducting 
currents. Only at extremely low temperatures, exponentially close to the absolute zero, the elec- 
trons could also form a gap leading to superconductivity of the whole system. In the present work 
we consider temperatures at which electrons are not condensed into Cooper pairs, and ignore the 
finite temperature effects. 

^Notc that for spin-dependent interactions the same effects of the charged condensate would 
give a generalization of the Ruderman-Kittel-Kasuya-Yosida (RKKY) potential [25] . 
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3. White Dwarf Cooling 



Most common dwarf stars have masses approximately equal to half of the mass 
of the Sun, central densities ~ (10^ — 10^) gjcw? and are composed mainly of carbon 
and/or oxygen. Above the largest of the three temperature scales - Tcryst, and 
Tc - the coohng process is determined by thermodynamics of the classical Bose gas 
of the ions. For lower temperatures however, the state of the star is determined by 
the relative magnitude of these scales. 

For present purposes, it is more convenient to rewrite the expressions for the 
temperature scales in terms of the mass density p measured in / crv? 

Teryst ^ 10^^^^ Od^A- 10^^ K, = 3.5 • lOVV^^ K , (11) 

where the baryon number of an ion was assumed to equal twice the number of 
protons (^4 = 2Z) and F ~ 180 was set. 

Helium White Dwarfs 

White dwarfs composed of helium constitute a smaller sub-class of dwarf stars (see, 
[9, 8] are references therein). They exhibit best conditions for the charged condensa- 
tion. Most of helium dwarfs arc believed to be formed in binary systems, where the 
removal of the envelope off the dwarf progenitor red giant by its binary companion 
happened before helium ignition, producing a remnant that evolves to a white dwarf 
with a helium core. Other astrophysical mechanisms for formation of isolated helium 
white dwarfs may also be possible [27]. In any event, helium dwarf masses range 
from ~ 0.5 Mq down to as low as (0.18 — 0.19) Mq, while their envelopes arc mainly 
composed of hydrogen. In this work we'll only consider He WDs whose hydrogen 
envelopes are thin, where no thermonuclear reactions are taking place. Such WDs 
cool by radiating off the stored in them heat. 

Following [7] , we will consider an over-simplified model of white dwarf cooling. 
Our treatment of what actually is an involved process, with significant uncertainties 
due to the envelope composition and opacity, should not be expected to give quanti- 
tatively precise predictions. Nevertheless, our approach is good enough, as we'll see, 
to capture the main difference that arises in cooling of dwarfs with the condensed 
cores. 

For definiteness, we consider cooling of a reference helium star of mass M = 
0.5 Mq with the atmospheric mass fractions of the hydrogen, and heavy elements 
(metallicity) respectively equal to 

X ~ 0.99, Zm (0.0002 - 0.002) . (12) 

The lower value of the metaUicity — 0.0002 is appropriate for the recently 
discovered 24 He WDs in NGC 6397 [9], but for completeness, we consider a wider 
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range for this parameter. Table 1 illustrates different physical characteristics of such 
a star with the mass density ~ 10® gjcrr? . 

The equation of state of a white dwarf interior is well modeled by the polytropes 
[7]. The average density of a non-relativistically degenerate star - which the helium 
dwarfs are an example of - may be about 5 times less than the central density. For 
the sake of simplicity, we will neglect the nontrivial density profile when dealing with 
thermodynamics of the dwarf stars and consider their average and central densities 
to be roughly equal. This won't lead to a significant error in the cooling analysis 
due to the fact that the variation of the mass density is not dramatic in the part of 
a star which encompasses most of its mass. 

As seen from Table 1, critical temperature for the reference star significantly ex- 
ceeds the crystallization temperature, while Debye temperature exceeds both. This 
star provides appropriate conditions for the formation of the charged condensate in 
its core. 

Table 1: Values of physical quantities for a reference helium white dwarf 



Physical quantity 


Numerical value 


Electron number density 
Mass density 
Separation between atomic nuclei 
Crystallization temperature 
Debye temperature 
Critical temperature 


(0.13 MeVf 
10® g/cw? 
10=^ fm 
3 • 10^ X 
A-W K 
10® K 



In the Appendix we give a brief summary of the derivation of the luminosity-core 
temperature relation for white dwarfs (for details see, e.g., [7]). The definition of 
luminosity of a cooling star, combined with equation (29) of the Appendix, gives 
the expression for the cooling rate: 

-cj— = L = CAmuT^-'. (13) 

Here denotes the specific heat per ion, m„ is the atomic mass unit and the 
constant C is inversely proportional to the atmospheric (envelope) metallicity of a 
star. In general, variations in atmospheric composition could change the cooling age 
significantly. Although the subject of atmospheric opacities is involved, the Kramer's 
approximation endowed with the atmospheric composition given in (12) is good 
enough for our purposes of comparing cooling rates of crystallized and condensed 
dwarfs. 

The expressions for specific heats of different components of dwarf cores in cor- 
responding regimes are given in Table 2. 

To quantify the effects of the charged condensation on the cooling rate, we con- 
sider the ratio of cooling times for two identical helium dwarfs with and without the 
charged condensate in their interiors. 
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Table 2: Specific heats of different components of a white dwarf core 



State 


Specific heat 


Classical Bose gas 
Quantum crystal 
Nonrelativistic Fermi Gas 

Relativistic Fermi Gas 


leyr-* ( T \3l. 
5 

2 Ep '^B 
(37r2)§ keTi 

3 ri/3 '^B 



Integrating (13), we find the following expression for the cooling time of a star 
in the Mestel regime 



k 

tUe — 



CAuiu 



-{Tp - TP) + Z-^{Tp - TP) 



(14) 



where Tj and Tq denote the final and initial core temperatures. The first term in the 
bracket on the right hand side corresponds to cooling due to classical gas of the ions 
and the second term corresponds to the contribution coming from the nonrelativistic 
Fermi sea. The latter is sub-dominant in the range of final temperatures we are 
interested in (the factor Z in front of this term is due to Z electrons per ion). Since 
Tf <^ Tq, the age of a dwarf star typically doesn't depend on the initial temperature. 
Neglecting the fermion contribution, we find time that is needed to cool down to 
critical temperature Tf — T^ 

Where an order of magnitude interval in (15) is due to the interval in the envelope 
metallicity composition given in (12). We also find the corresponding luminosities 



M fT 



7/2 

L{T,) ~ (10« erg/s)^ (^^ j ~ 1.5 • (10"^ - 10-^)Lq , (16) 

which arc in the range of observable luminosities (L0 ~ 3.84 ■ 10'^'^ erg/ s). 

After the condensation, specific heat of the system dramatically drops as the 
collective excitations of the condensed nuclei become massive and "get extinct" . A 
contribution from the Fermi sea, which is strongly suppressed by the value of Fermi 
energy, becomes the dominant one. The phase transition itself would take some time 
to complete, and the drop-off in specific heat will not be instantaneous. During that 
time a heat-transfer from the bosonic sector to the fermionic one will take place, but 
this will only change temperature of the fermions by a factor of (1 -|- 1/Z), which, 
given our approximations, can be ignored. Effects of finite duration of the phase 
transition on the shape of the luminosity function will be discussed in the next 
section. Here, for simplicity we approximate the transition to be instantaneous, and 
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retain only the fermion contribution to specific heat below Tc. Then, the expression 
for the age of the star for Tf < T^, reads as follows 



t 



CArriu 



k 



B 



3 Ef 



To ' . 



(17) 



Notice the difference of (17) from (14) - in the former Tf < Tc and it is Tj that 
enters as final temperature in the fermionic part, while Tc should be taken as the 
final temperature in the bosonic part. 

20 



15 



^10 



400000 



600000 800000 
Tf (K) 



1000000 



Figure 1: The ratio of the ages, as a function of final temperature, for two identical 
helium dwarf stars, with and without the interior condensation. 

Figure 1 gives the ratio of ages, t] = tHe/i'nei identical helium dwarf 

stars, with and without the interior condensation. The considered range of final 
temperature is (3 ■ 10^ - 10^) K. At the upper edge of this interval, the luminosity 
of the star, depending on the metallicity interval given in (12), is (10~^ — 1O~^)L0, 
while for the lower edge of the temperature interval we get 2 ■ (10~^ — 10~^)Lq. As 
is seen from this figure, the charged condensation substantially increases the rate of 
cooling. At the lower end of the considered temperature interval the age is about 
twenty times less than it would have been without the condensation phase. 

2 

With the decrease of density, the critical and Debye temperatures drop as oc pa 
and oc p2 respectively, whereas the crystallization temperature drops only as 
Tcryst OC p3. Heuce, the critical and crystallization temperatures would become 
equal for the dwarfs of low density. This would happen for densities ~ 2 ■ 10^ g/cm?, 
which are very low and unlikely to be present in white dwarfs cores. Hence, we'd 
expect that all He WDs should have a large fractions of their cores in the condensed 
state. 



13 



To summarize, the charged condensation of white dwarf interiors would have a 
significant impact on coohng curves of hehum dwarf stars at the faint end of the 
luminosity spectrum. Starting from a certain temperature scale, ~ Tc , the ages 
of the stars become effectively " frozen" , whereas the temperature and luminosity 
decrease rather fast. Such a change in the cooling regime would leave its signatures 
on the luminosity function of helium white dwarfs, which we'll discuss below. 



Luminosity Function 

The himinosity function (LF) (jyijr) '^^'^ nsed to test the white dwarf evolution 
models against the observational cfata; it is defined as follows: 

(f){^)d{log{-^)) = number-density of WD's per unit interval of (log(^)). (18) 

Under the assumption of the standard star formation rate, uniform both in space 
and time, the luminosity function takes a simple power-law form 



(j) oc 



dt 



n(k + l) , 

ocL 4 -1, (19) 



where the constants k and n are the exponents defining temperature dependence of 
specific heat and luminosity: 

c^oc^^ LocTI (20) 

For instance, n equals to 8/7, once the Kramer's opacity is adopted for the descrip- 
tion of the dwarf atmosphere. The log of the LF corresponding to cooling in the 
Mestel regime (c^ = 3/^^/2), is just a line with a slope equal to 5/7. 

The condensation of the core would induce significant deviations from the Mestel 
curve for helium white dwarfs. The exponential suppression of specific heat of the 
ions in the condensate doesn't occur instantaneously. Typically, for a quantum 
liquid, there is a transition region of width ~ ATp, in which the specific heat changes 
in a certain way. Therefore, starting from some temperature Tc — ATc, at which 
the specific heat of the ions can be neglected and a stage of rapid cooling begins 
(the age of the star "freezes"), the luminosity function starts to drop dramatically 
(Fig. 2). The stars, previously belonging to the same logarithmic luminosity bin, 
will spread among a few bins as a result of suddenly accelerated cooling. 

Not much can be said with certainty of ATc and the temperature dependence 
of specific heat near the phase transition point. In a well-known case of quantum 
liquid of neutral helium atoms for instance, the specific heat exhibits the so-called 
"Lambda" behavior. In analogy with the latter, one may expect the specific head to 
increase near the phase transition point before dropping dramatically below T(.. If so, 
the parameter ATc would give the width of the transitional region. Depending on the 
details of the function c^(T) near the phase transition point, the luminosity function 
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'Sun 



Figure 2: A schematic sketch of the luminosity function for hehum-core white dwarfs. 
The absolute normalization of the function is set by the constant b which is deter- 
mined by their formation rate. The blue line represents the Mestel regime. The 
shape of the luminosity function near the condensation point depends on the details 
of the corresponding phase transition. 
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of He WDs could take different shapes near its maximum, as it is schematically 
illustrated in Fig. 2. 

What is independent of the above uncertainties, however, is the fact that the LF 
will experience a significant drop-off after the charged condensation phase transition 
is complete. This is due to the "extinction" of the bosonic quasiparticles below the 
phase transition point. In fact, the LF function will drop by a factor of ~ 200, 
which would be reflected as a drop-off by about two units of the vertical axis from 
its maximum shown on Fig 2. After the drop-off the log of the LF function starts 
growing again but now with the slope equal to 3/7 - due to the specific heat of the 
electrons. The latter part of the LF is shown only for one of the curves on Fig 2. 
Whether this drop-off can explain the termination of the 24 He WD sequence found 
in [9], remains to be seen in more detailed studies. 



Superdense Carbon White Dwarfs 

Finally, we briefly mention another possible subclass of white dwarf stars, that could 
undergo the core condensation. These are superdense dwarfs composed mainly of 
carbon, with masses close to the Chandrasekhar limit ~ 1.4 Mq, and central densi- 
ties ~ 10^*^ 5'/cm'^. Such a star would be very close to the neutronization threshold 
(5- lO^^g/cm^ for carbon nuclear matter), and it's critical temperature ~ 8-10'^ K 
would be greater than the crystallization temperature T^ryst ~ 4-10^ K, making con- 
densation a possibility^. 

Assuming the helium dominated envelope {Y denotes the helium fraction) with 
a 10% metallicity 

a:~o, y~o.9, z„i-o.i, (21) 

we find the cooling time needed to reach down to temperature Tc ~ 8 • 10^ K 

This temperature corresponds to luminosity L ~ 3.4Lq. Once the star enters the 
condensation regime, it cools faster. For instance, it would reach the luminosity 
L ~ 10-U in t = 0.7 Gyr after the condensation. On the other hand, if one 
assumes that this star crystallizes instead of condensing its core, then it would reach 
the crystallization point in 10~^ Gyr, while it would take only ~ 0.1 Gyr longer for 
this star to cool down to L ~ 1O~^L0, as compared with the condensation case. 
In the crystallized phase for this star the cooling due to the fermion contributions 
is significant. Therefore, the log of the LF for such a star would not have a sharp 



''At these high temperatures the WD coohng rate is significantly affected by the neutrino emis- 
sion. Our goal, however, is to identify qualitative differences of cooling in the condensate and 
crystalline phases. For this purpose, and for simplicity, we ignore the neutrino contributions, 
which should certainly be taken into account for these WDs in more precise studies. 
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drop-off; instead, after the phase transition the LF would just change its slope to 
3/7 in the low luminosity region. 

We also note that for dwarf stars with relativistically degenerate cores, the poly- 
trope models yield more dramatic density profiles - density in the center exceeds the 
average density by a factor of ~ 50. Since the crystallization/critical temperature 
equality is achieved at densities ~ 10^ g/cm? for carbon interiors, we should not 
expect the entire core to condense. In realistic calculations of cooling rates of the 
superdense carbon WDs this should also be taken into account. 

Magnetic properties of condensed WDs, and their similarity to type II supecon- 
ductors, will be discussed in [28]. 

Finally, we mention that charged condensation may also take place in some other 
astrophysical objects where the densities and temperatures are appropriate. Crusts 
of neutron stars may be a place to look at. 
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Appendix 

For purposes of studying their cooling, white dwarfs are well described by a simple 
2-component model. The core consists of the ion gas and degenerate electrons with 
large mean free path and heat conductivity, making it almost isothermal. The 
coohng takes place through the non-degenerate surface layer, the envelope, that 
surrounds the core. The photon diffusion equation, describing the energy fiow from 
the core to the outer layers, has the following form 

where L is luminosity, aT^- energy density of a blackbody, n is the opacity of the 
stellar matter and p-mass density. We can rewrite the last equation in the following 
form 

^ ^ ^'^Ji_ C24^ 

dr AacT^A-nr'^' ^ ' 

For the opacity, we use Kramer's approximation k = KopT~^-^. 
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The equation of the hydrostatic equilibrium (where P denotes the pressure and 
m(r) - the mass inside a sphere of radius r) 



dP ^ _Grnpp^ 
dr 

combined with the equation of state for a classical ideal gas 

P = -^kBT , (26) 

(where ^ and rriu are the mean molecular weight and atomic mass unit respectively) , 
can be integrated with the boundary condition P = at T = to obtain 



2 AacA-nGM jjirnu 
B.5 3 KqL ks 



2 



(27) 



where Kq = 4.34 • 10^'^ ^^(1 + X) cw? / g is the commonly used value for the opacity 
constant, X is the mass fraction of hydrogen, and Z.^ - that of heavy elements (all 
elements except hydrogen and helium), while M is the mass of the star. Equation 
(27) works as long as we deal with nondegenerate matter. To estimate the limits of 
its applicability, we should equate the pressure of nondegenerate electrons (obtained 
from (26) by replacing by mean molecular weight per electron //g) to that of the 
nonrelativistic electron gas in the outer layers of the core 

P*kBT^ _ 3t7rt Ifp'J^ , . 

where and T* are the density and temperature of the star at the core-surface 
boundary. Combining equations (28) and (27), we obtain a useful relation for the de- 
pendence of the luminosity on the temperature at the core-surface boundary (which, 
due to the isothermahty of the core, can be used as a good characteristic of the core 
temperature) 

lli Zm{l + X) Mq 

where 7 ~ 5.7 ■ 10^ erg/ s is a universal constant, and C is a constant that depends 
on the chemical composition of the envelope and varies from a star to star. A more 
general form of the last equation may be written as follows: 

LocTn, (30) 

where n is a constant, equal to 8/7 in case of the envelope with Kramer's opacity. 
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